/*
 #
 #  Files       : SmallMatrixIO.h
 #                ( C++ header file )
 #
 #  Description : The SmallMatrix Library
 #                ( http://code.google.com/p/smallmatrix )
 #
 #  Copyright   : Olivier Juan
 #                ( http://www.mas.ecp.fr/vision/Personnel/juan/ )
 #
 #  License     : CeCILL-C
 #                ( http://www.cecill.info/licences/Licence_CeCILL-C_V1-en.html )
 #
 #  This software is governed by the CeCILL-C license under French law and
 #  abiding by the rules of distribution of free software.  You can  use,
 #  modify and or redistribute the software under the terms of the CeCILL-C
 #  license as circulated by CEA, CNRS and INRIA at the following URL
 #  "http://www.cecill.info".
 #
 #  As a counterpart to the access to the source code and  rights to copy,
 #  modify and redistribute granted by the license, users are provided only
 #  with a limited warranty  and the software's author,  the holder of the
 #  economic rights,  and the successive licensors  have only  limited
 #  liability.
 #
 #  In this respect, the user's attention is drawn to the risks associated
 #  with loading,  using,  modifying and/or developing or reproducing the
 #  software by the user in light of its specific status of free software,
 #  that may mean  that it is complicated to manipulate,  and  that  also
 #  therefore means  that it is reserved for developers  and  experienced
 #  professionals having in-depth computer knowledge. Users are therefore
 #  encouraged to load and test the software's suitability as regards their
 #  requirements in conditions enabling the security of their systems and/or
 #  data to be ensured and,  more generally, to use and operate it in the
 #  same conditions as regards security.
 #
 #  The fact that you are presently reading this means that you have had
 #  knowledge of the CeCILL-C license and that you accept its terms.
 #
*/

#ifndef _SMALLMATRIXIO_H
#define _SMALLMATRIXIO_H

#include <SmallMatrix/SmallMatrix.h>
#include <iostream>
#include <sstream>

#pragma warning (disable : 4127)
template <typename T> INLINE std::ostream& operator<<(std::ostream& out,const SmallMatrix::Matrix_Base<T>& vect) {
	if (T::COL == 1) {
		out << "(";
		for (int i=0;i<T::ROW;i++) {
			out<<vect(i,0);
			if (i!=(T::ROW-1)) out<<", ";
			else out<<")"<<std::endl;
		}
	} else {
		for (int i=0;i<T::ROW;i++) {
			for (int j=0;j<T::COL;j++) {
				out<<vect(i,j);
				if (j!=(T::COL-1)) out<<"\t";
			}
			out<<std::endl;
		}
	}
	return out;
}

template <typename T> INLINE std::istream& operator>>(std::istream& in, SmallMatrix::Matrix_Base<T>& vect) {
	for (int i=0;i<T::SIZE;i++) {
		in>>vect[i];
	}
	return in;
}

#pragma warning ( disable : 4100)
template <typename T> std::string AfficheType(const SmallMatrix::Matrix_Base<T> & A) {
	T temp;
	return AfficheType(temp);
}

template <int M,int N> std::string AfficheType(const SmallMatrix::Matrix<M,N> & A) {
	std::stringstream sstr;
	sstr<<"Matrix "<<M<<"x"<<N;
	return sstr.str();
}

template <int N> std::string AfficheType(const SmallMatrix::SquareMatrix<N> & A) {
	std::stringstream sstr;
	sstr<<"SquareMatrix "<<N<<"x"<<N;
	return sstr.str();
}

template <int N> std::string AfficheType(const SmallMatrix::SymMatrix<N> & A) {
	std::stringstream sstr;
	sstr<<"SymMatrix "<<N<<"x"<<N;
	return sstr.str();
}

template <int N> std::string AfficheType(const SmallMatrix::Vecteur<N> & A) {
	std::stringstream sstr;
	sstr<<"Vecteur "<<N;
	return sstr.str();
}

template <typename Functor, typename T, typename TT> INLINE std::ostream& operator<<(std::ostream& out,const SmallMatrix::Expressions<Functor, T, TT>& Expr) {
	typedef typename SmallMatrix::Expressions<Functor, T, TT>::result_type result_type;
	if (result_type::COL == 1) {
		out << "(";
		for (int i=0;i<result_type::ROW;i++) {
			out<<Expr(i,0);
			if (i!=(result_type::ROW-1)) out<<", ";
			else out<<")"<<std::endl;
		}
	} else {
		for (int i=0;i<result_type::ROW;i++) {
			for (int j=0;j<result_type::COL;j++) {
				out<<Expr(i,j);
				if (j!=(result_type::COL-1)) out<<"\t";
			}
			out<<std::endl;
		}
	}
	return out;
}

#pragma warning ( default : 4127)

#endif
